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Abstract 

We compute the corrections to the orbital Lense-Thirring effect (or frame-dragging) 
in the presence of spacetime torsion. We analyze the motion of a test body in the 
gravitational field of a rotating axisymmctric massive body, using the parametrized 
framework of Mao, Tegmark, Guth and Cabi. In the cases of autoparallel and extremal 
trajectories, we derive the specific approximate expression of the corresponding system 
of ordinary differential equations, which are then solved with methods of Celestial Me- 
chanics. We calculate the secular variations of the longitudes of the node and of the 
pericenter. We also show how the LAser GEOdynamics Satellites (LAGEOS) can be 
used to constrain torsion parameters. We report the experimental constraints obtained 
using both the nodes and perigee measurements of the orbital Lense-Thirring effect. 
This makes LAGEOS and Gravity Probe B (GPB) complementary frame-dragging and 
torsion experiments, since they constrain three different combinations of torsion param- 
eters. 

Keywords: Riemann-Cartan spacetime, torsion, autoparallel trajectories, frame dragging, 
geodetic precession, satellite laser ranging, Gravity Probe B. 

1 Introduction 

In recent years a lot of effort has been devoted to measure gravitomagnetic effects due to 
Earth's rotation [I], [2], [3] predicted by the theory of General Relativity (GR). In particular, 



the Lense-Thirring effect on the orbital motion of a test body can be measured by using 
the satellite laser ranging (SLR) technique, whose data are provided by the ILRS^j By 
analyzing the laser ranging data of the orbits of the satellites LAGEOS and LAGEOS II, a 
measurement of the Lense-Thirring effect was obtained by Ciufolini and Pavlis [JJ. 
SLR missions can also be useful to test modifications of GR, such as torsion theories of 
gravity. A class of theories allowing the presence of torsion is based on Riemann-Cartan 
spacetime, which is endowed with a metric g^ u and a compatible connection. The resulting 
connection v turns out to be nonsymmetric, and therefore it originates a non-vanishing 
torsion tensor. We refer to [S], [Hj for the details. 

In standard torsion theories the source of torsion is considered to be the intrinsic spin of 
matter [5], [6], [7], [8], which is negligible when averaged over a macroscopic body. Therefore 
spacetime torsion would be observationally negligible in the solar system. Nevertheless, in 
j9] Mao, Tegmark, Guth and Cabi (MTGC) argue that the presence of detectable torsion 
in the solar system should be tested experimentally, rather than derived by means of a 
specific torsion model. For this reason, in [9] a theory-independent framework based on 
symmetry arguments is developed, and it is determined by a set t±, t2, w%, . . . , W5 of seven 
parameters describing torsion and three further parameters J-,Q,T~L describing the metric. 
Here, by theory-independent framework, we mean the following: the metric and the con- 
nection are parametrized, around a massive body, with the help of symmetry arguments, 
without reference to a torsion model based on a specific Lagrangian (or even on specific 
field equations). 

This parametrized framework can be used to constrain t\ , ti , W\ , . . . , W5 from solar system 
experiments. In particular, MTGC suggest that GPB [10J is an appropriate experiment for 
this task, and in [9] they compute precessions of gyroscopes and put constraints on torsion 
parameters from GPB measurements. In [11] Hehl and Obukhov argue that measuring 
torsion requires intrinsic spin, and criticize the approach of MTGC, since GPB gyroscopes do 
not carry uncompensated elementary particle spin. Nevertheless, we accept the general idea 
that the precise form of the coupling of torsion to matter should be tested experimentally, 
and that actual experimental knowledge leaves room for nonstandard torsion theories which 
could yield detectable torsion signals in the solar system. In the present paper we apply 
the parametrized framework developed by MTGC for the computation of satellites orbits 
around Earth and we put a different set of constraints on torsion parameters from SLR 
measurements. 

MTGC also address the question of whether there exists a specific gravitational Lagrangian 
fitting in the parametrized framework and yielding a torsion signal detectable by the GPB 
experiment. As an example they quote the theory of Hayashi and Shirafuji (HS) in |12j 
where a massive body generates a torsion field, and they propose what they call the Einstein- 
Hayashi-Shirafuji (EHS) Lagrangian, interpolating GR and HS Lagrangians in a linear way. 
However, MTGC consider only a gravitational Lagrangian in vacuum, so that they cannot 
derive the equations of motion of test bodies from the gravitational field equations, which 
would require a suitable matter coupling. 

The EHS model has been criticized by various authors. In the paper [T3], Flanagan and 
Rosenthal show that the linearized EHS theory becomes consistent only if the coefficients 
in the Lagrangian are constrained in such a way that the resulting predictions coincide with 
those of GR. In the paper [TJJ , Puetzfeld and Obukhov derive the equations of motion in the 
framework of metric- affine gravity theories, which includes the HS theory, and show that 
only test bodies with microstructure (such as spin) can couple to torsion. In conclusion, 
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the EHS theory does not yield a torsion signal detectable for GPB. For these reasons, in j9] 
the EHS Lagrangian is proposed not as a viable physical model, but as a pedagogical toy 
model fitting in the parametrized framework, and giving an illustration of the constraints 
that can be imposed on torsion by the GPB experiment. In the present paper we will not 
consider such a toy model. 

As also remarked by Flanagan and Rosenthal in |13| . the failure of constructing the spe- 
cific EHS Lagrangian does not rule out the possibility that there may exist other torsion 
theories which could be usefully constrained by solar system experiments. Such torsion 
models should fit in the above mentioned theory-independent framework, similarly to a 
parametrized post-Newtonian framework including torsion. We remark that the parametrized 
formalism of MTGC does not take into account the intrinsic spin of matter as a possible 
source of torsion, and in this sense it cannot be a general torsion framework. However, it 
is adequate for the description of torsion around macroscopic massive bodies in the solar 
system, like planets, being the intrinsic spin negligible when averaged over such bodies. 
Therefore we think it is worthwhile to continue the investigation of observable effects in 
the solar system of nonstandard torsion models within the MTGC parametrized formalism, 
under suitable working assumptions. In particular, our aim is to extend the GPB gyroscopes 
computations made in [9] to the case of motion of satellites. 

In the present paper we compute the corrections to the orbital Lense-Thirring effect due to 
the presence of spacetime torsion described by ti, £2, Wi, . . . , w§. We consider the motion of 
a test body in the gravitational field of a rotating axisymmetric massive body, under the 
assumption of slow motion of the test body. Since we use a parametrized framework without 
specifying the coupling of torsion to matter, we cannot derive the equations of motion of 
test bodies from the gravitational field equations. Therefore, in order to compute effects 
of torsion on the orbits of satellites, we will work out the implications of the assumption 
that the trajectory of a test body is either an extremal or an autoparallel curve. Such 
trajectories do not need to coincide when torsion is present. 

As in the original paper of Lense and Thirring [15] , we characterize the motion using the six 
orbital elements of the osculating ellipse. In terms of these orbital elements, the equations 
of motion then reduce to the Lagrangian planetary equations. We calculate the secular 
variations of the longitude $7 of the node and of the longitude uj of the pericenter. The 
computed secular variations show how the corrections to the orbital Lense-Thirring effect 
depend on the torsion parameters, and it turns out that the dependence is only through 
w%, . . . , W5. The data from the LAGEOS satellites are then used to constrain the relevant 
linear combinations of the torsion parameters. More precisely, we constrain two different 
linear combinations of wi, . . . , w$ by using first the measurements of the nodes of LAGEOS 
and LAGEOS II, and then the measurements of the nodes of LAGEOS and LAGEOS II 
and of the perigee of LAGEOS II. In particular, torsion parameters cannot be constrained 
by satellite experiments in the case of extremal trajectories. 

While the torsion perturbations to the Lense-Thirring effect depend only on wi, . . . , w$, 
it turns out that another relevant relativistic effect, namely the geodetic precession (or de 
Sitter effect), depends on the parameters t\ and t2, and on a further parameter £3. This 
latter parameter is involved in a higher order parametrization of torsion, which is necessary 
for the description of the geodetic precession effect, while it is not necessary at the order of 
accuracy required in the present paper. All computations of orbital geodetic precession with 
torsion of a satellite are performed in the companion paper |16] , to which we will sometimes 
refer for details. 

The paper is organized as follows. In Section [2] we briefly recall the notion of spacetime with 
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torsion. In Section[3]we discuss the case of extremal trajectories. In Section|4]we analyze the 
equations of autoparallel trajectories and derive the related system of ordinary differential 
equations to first order. The expression of the system clearly reveals the perturbation due 
to torsion with respect to the Lense-Thirring equations. In Section [5] we derive the time 
evolution of the orbital elements, by applying the classical perturbation theory of Celestial 
Mechanics, in particular the Gauss form of the Lagrange planetary equations. In Section 
[6] we calculate the secular variations of the orbital elements. In Section [7] we recall some 
results from [16] where torsion solar perturbations are computed. These results will be used 
in Section [8j where we give the observational constraints that the LAGEOS experiment can 
place on torsion parameters. Conclusions are drawn in Section [9| For convenience of the 



reader, in the appendix (Section 10) we recall from [9] how to parametrize the metric and 
torsion tensors, and hence how to parametrize the connection, under suitable symmetry 
assumptions. 

2 Spacetime with torsion 

A manifold equipped with a Lorentzian metric g^ v and a connection V\ v compatible with 
the metric is called a Riemann-Cartan spacetime [5], |BJ. Compatibility means that V ^g v \ = 
0, where V denotes the covariant derivative. We recall in particular that for any vector field 

The connection is determined uniquely by g^ u and by the torsion tensor 

c X / "pA -pA ^ 

as follows: 

= (21 » 

where {•} is the Levi-Civita connection, defined by 

/it} = ^ 9up + Qv9w ~ dp9 ^ ' (2 ' 2) 
and 

■K-nv = ~S fJLV — S Ufi — S M „ (2-3) 

is the contortion tensor. In the particular case when T \ v is symmetric with respect to 
fi, v the torsion tensor vanishes. We will be concerned here with the case of nonsymmetric 
connections T^„ v . The case of vanishing torsion tensor corresponds to Riemann spacetime of 
GR, while the case of vanishing Riemann tensor corresponds to the Weitzenbock spacetime 



In the present paper we use the natural gravitational units c = 1 and G = 1. We will 
assume that Earth can be approximated as a uniformly rotating spherical object of mass m 
and angular momentum J. Following [9], we use spherical coordinates (r, 9, (f>) for a satellite 
moving in the gravitational field of Earth, and we introduce the dimensionless parameters 
e m = m/r and ej = J / (mr). Since the radii of the LAGEOS orbits (about 6000 km altitude) 
are much larger than Earth's Schwarzschild radius, it follows that e m « 1. Moreover, since 
Earth is slowly rotating, we have ej << 1. Therefore, all computations will be carried out 
perturbatively to first order in e m and ej. 
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Under spherical axisymmetry assumptions, the metric tensor g^ u and the torsion tensor 
p have been parametrized to first order in [9] . Accordingly, is parametrized by three 
parameters H, Q, and S fll , p is parametrized by seven parameters t%, t2, wi, . . . , W5, 

s ^u X = S fiu X (h,t2,wi, • • .,w 5 ,r,6,</>) . 

Therefore r A ^„ becomes an explicit function of all metric and torsion parameters. It turns 
out that ti,t% contribute to geodetic precession, while w\ , • • • , 105 contribute to the frame- 
dragging precession. In the Appendix we report the explicit expressions of the parametrized 
metric and torsion tensors, and of the connection, that will be needed in the sequel of the 
paper. 



3 Equations of extremal trajectories 

In GR structureless test bodies move along geodesies. In a Riemann-Cartan spacetime 
there are two different classes of curves, autoparallel and extremal curves, respectively, 
which reduce to the geodesies of Riemann spacetime when torsion is zero [5] . Autoparallels 
are curves along which the velocity vector is transported parallel to itself by the connection 
r A MJy . Extremals are curves of extremal length with respect to the metric g^ u . The velocity 
vector is transported parallel to itself along extremal curves by the Levi-Civita connection. 
In GR the two types of trajectories coincide while, in general, they may differ in presence of 
torsion. They are identical when the torsion is totally antisymmetric [5], a condition which 
is not satisfied within our parametrization. 

The equations of motion of bodies in the gravitational field follow from the field equations 
due to the Bianchi identities. The method of Papapetrou (TTj can be used to derive the 
equations of motion of a test body with internal structure, such as for instance a small 
extended object that may have either rotational angular momentum or net spin. In standard 
torsion theories the trajectories of test bodies with internal structure, in general, are neither 
autoparallels nor extremals [5], [S], [IB]) while structureless test bodies, such as spinless test 
particles, follow extremal trajectories. 

The precise form of the equations of motion of bodies in the gravitational field depends 
on the way the matter couples to the metric and the torsion in the Lagrangian (or in 
the gravitational field equations). As explained in the Introduction, we do not specify a 
coupling of torsion to matter, hence we do not specify the field equations. Moreover, in 
our computations of orbits of a satellite (considered as a test body), we will neglect its 
internal structure. In a theory-independent framework we cannot derive the equations of 
motion from the gravitational field equations, hence we need some working assumptions 
on the trajectories of structureless test bodies: we will investigate the consequences of 
the assumption that the trajectories are either extremal or autoparallel curves. Assuming 
the trajectory to be an extremal is natural and consistent with standard torsion theories. 
However, extremals depend only on the parameters of the metric, so that new predictions 
related to torsion cannot arise. We will quickly report the computations for the sake of 
completeness, since the metric parameters can be immediately related to the Parametrized 



Post Newtonian (PPN) parameters (see (10.2)), and the orbital Lense-Thirring effect in the 
case of extremal trajectories and a PPN metric is known. 
The system of equations of extremal trajectories reads as 

d 2 x x f A 1 dx^ dx u , , 
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where r is the proper time. For slow motion of the satellite we can make the substitution 
dr ~ dt, so that 

d 2 x a { a "I dx^ dx u 
jiv) dt dt 



dt 2 



+ 



0. 



for a £ {1,2,3}. We assume that the velocity of the satellite is small enough so that we 
can neglect the quadratic terms in the velocity. Then, being x° = t we have 



d 2 x° 



+ 



+ 2 



a "1 dx? 
0/3 J lit 



0, 



(3.2) 



for (3 £ {1,2,3}. 

All perturbations considered here are so small that can be superposed linearly. Since we 
are only interested in the perturbations due to Earth's rotation, as in the original Lense- 
Thirring paper |15j we are allowed to neglect the quadratic terms in the velocities which 
yield an advance of the perigee of the satellite. The value of the advance of the perigee for 
an extremal orbit and a PPN metric can be found in [2, Chapter 7, formula (7.54)]. 
We use for x a spherical coordinates (r, 9, 4>). The Levi-Civita connection {•} can be obtained 
from the expression of T „ v given in the Appendix by setting to zero all torsion parameters 



ti, t%, W\, . . . , W5. Substituting the resulting expression in (3.2 ) one gets 



rr - (H/2)e m + Q(f>r sin 2 9e m ej = 0, 



2<5</> sin # cos 6>e m e j = 0, 



(3.3) 



hr 2 sin 



— Qr sin 0e m ej + IQOr cos 9e m ej = 0. 



The equations of motions (3.3 ) depend neither on the metric parameter T nor on the torsion 



parameters. System (3.3) to lowest order becomes 

dv T~L m „ 
dt = 2"^ 6r ' 

where e r is the unit vector in the radial direction. Imposing the Newtonian limit yields 
T~L = —2 as in a PPN metric (see also ^ formula (23)]). 
We now transform (3.3) in rectangular coordinates x = rsinl 
z = rcos( 



1 cos < 



y 



r sin v sin < 



We compute the second derivatives of x,y,z with respect to time in the 
approximation of slow motion. Neglecting all terms containing squares and products of 
first derivatives with respect to (r, 8, </)), we get 



f sin 9 cos d> + 9r cos 9 cos 1 



or sin t) sin 1 



y = r sin 9 sin <p + 9r cos 9 sin <j) + <pr sin 9 cos 4>, 



, z = r cos 9 — 9r sin 9. 



(3.4) 



Using (3.3) and (3.4) we obtain the following system for the equations of motion: 



r 

Cm 



n 

-x — y 



y + Q 



[x 2 + y 2 - 2z 2 ) y + 3yzz 
(x 2 + y 2 — 2z 2 ) x + 3xzz 



(3.5) 



z + G^^-3z (yx - xy) . 
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Note that when Q = —2 system (3.5) reduces to the equations of motion found by the Lense- 
Thirring [15] formula (15)]. Hence the relativistic perturbation of the Newtonian force is 
just multiplied by the factor —Q/2 with respect to the original Lense-Thirring equations. It 
follows that the formulae of precession of the orbital elements of a satellite can be obtained 
by multiplying the original Lense-Thirring formulae |15|. formula (17)] by the factor —Q/2. 
The details of the computation, based on the Lagrange planetary equations of Celestial 
Mechanics, can be also retrieved from the computations for autoparallel trajectories given 
in the next sections, by setting to zero all torsion parameters t\, t2, Wi, . . . , w§. 
Using the standard astronomical notation, we denote by the longitude of the node and 
by u) the argument of the perigee of the satellite's orbit. The secular contributions to the 
variations of f2 and u are: 

( <5fi ) sec = ~ a 3(l_ e 2)3/2 *' ^ scc = fl 3 (1 _ e 2)3/2 *> ( 3 " 6 ) 

where a is the semimajor axis of the satellite's orbit, e is the eccentricity, i is the orbital 



inclination, and t is time. When Q = —2 the quantities in (3.6) reduce to the classical 
corresponding Lense-Thirring ones. 

Since the expressions of (<5J7) S0C and (5u>) sec depend only on Q, the measurements of satellites 
experiments cannot be used to constrain the torsion parameters. 



4 Equations of autoparallel trajectories 

In standard torsion theories the trajectories of structureless test bodies follow extremal 
trajectories [5], [6], which depend only on the metric. However, new predictions related to 
torsion may arise when considering the autoparallel trajectories. In the following we give 
some motivations which make worthwhile the investigation of autoparallel trajectories. 
Since in spacetime with torsion parallelograms are in general not closed, but exhibit a 
closure failure proportional to the torsion, Kleinert and Pelster argue in |19j that the vari- 
ational procedure in the action principle for the motion of structureless test bodies must 
be modified. In the standard variational procedure for finding the extrema of the action, 
paths are varied keeping the endpoints fixed in such a way that variations form closed paths. 
However, in the formalism of [19] , the closure failure makes the variation at the final point 
nonzero, and this gives rise to a force due to torsion. When this argument is applied to 
the action principle for structureless test bodies it turns out that the resulting torsion force 
changes extremal trajectories to autoparallel ones (see |19| for the details). Kleinert and 
Shabanov find an analogous result in [20] where they show that the geometry of spacetime 
with torsion can be induced by embedding its curves in a euclidean space without torsion. 
Kleinert et al. also argue in |19j . j20j that autoparallel trajectories are consistent with the 
principle of inertia, since a structureless test body will change its direction in a minimal 
way at each time, so that the trajectory is as straight as possible. 

The approach of Kleinert et al. has been criticized by Hehl and Obukhov in [11] since the 
equations of autoparallel trajectories have not been derived from the energy- momentum 
conservation laws. Kleinert investigates this issue in |21] and finds that, due to the closure 
failure, the energy-momentum tensor of spinless point particles satisfies a different conser- 
vation law with respect to the one satisfied in torsion theories such as [5] , [6] . The resulting 
conservation law yields autoparallel trajectories for spinless test particles. Kleinert then ad- 
dresses the question of whether this new conservation law allows for the construction of an 
extension of Einstein field equations to spacetime with torsion. The author gives an answer 
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for the case of torsion derived from a scalar potential (see [6] for a discussion of this kind 
of torsion). In this case the autoparallel trajectories are derived from the gravitational field 
equations via the Bianchi identities, though the field equation for the scalar field, which is 
the potential of torsion, is unknown. 

In [22] Dereli and Tucker show that the theory of Brans-Dicke can be reformulated as a field 
theory on a spacetime with dynamic torsion determined by the gradient of the Brans-Dicke 
scalar field. Then in [23] they suggest that the autoparallel trajectory of a spinless test 
particle in such a torsion geometry is a possibility that has to be taken into account. In |23j 
the autoparallel trajectories of massive spinless test particles are analyzed in the background 
of a spherically symmetric, static solution to the Brans-Dicke theory and the results are 
applied to the computations of the orbit of Mercury. In [23] the autoparallel trajectories 
of spinless particles are analyzed in the background of a Kerr Brans-Dicke geometry. In 
[25], [26] the equations of autoparallel trajectories are derived from the gravitational field 
equations and Bianchi identities, in the special case of matter modeled as a pressureless 
fluid, and torsion expressed solely in terms of the gradient of the Brans-Dicke scalar field. 
The above quoted results show that there is an interest in the autoparallels in spacetime 
with torsion, which make worthwhile their investigation in the present paper. The system 
of equations of autoparallels reads as 

dr 2 + dr dr ' 1 ' 

where r is the proper time [27]. Observe that only the symmetric part ^(r^ + T A ) of 



the connection enters in (4.1); moreover, starting from (4.1) the totally antisymmetric part 
of S\^ u cannot be measured. 

The trajectory of a test body has to be a time-like curve. Since the connection is compatible 
with the metric, the quantity 9fiu^r^§f is conserved by parallel transport. The tangent 
vector to the trajectory undergoes parallel transport by the connection along the au- 
toparallel. Therefore, an autoparallel that is time-like at one point has this same orientation 
everywhere, so that the trajectory is strictly contained in the light cone determined by g^, 
in a neighbourhood of every of its points. Hence the compatibility of the connection with 
the metric ensures that autoparallels fulfil a necessary requirement for causality. 
For slow motion of the satellite we can make the substitution dr ~ dt, so that 

d 2 x a dx^ dx v 

~df r+ ^ dt dt ' 

for a £ {1, 2, 3}. Again, we assume that the velocity of the satellite is small enough so that 
we can neglect the terms which are quadratic in the velocity. Then, being x° = t we have 

+ r° 00 + (r> + r»— = o, (4.2) 

for (3 £ {1,2,3}. 

As in the previous section, all the perturbations that we are considering here are so small 
that can be superposed linearly. We are allowed to neglect the quadratic terms in the 
velocities which yield an advance of the perigee of the satellite. Such an advance of the 
perigee for an autoparallel orbit in presence of torsion has been computed in |16J . 



We use for x a spherical coordinates (r,8,(j)). Substituting in (4.2) the expression of T 



fJ,U 
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given in the Appendix one gets 

rr + Ce m + Vtfir sin 2 9e m ej = 0, 

Or — Bcfi sin 6 cos 9e m ej = 0, 

4>r 2 sin 6 + «4r sin 0e m ej + ,6#r cos 0e m ej = 0, 



(4.3) 



where 



A = -Q + w\ - W3, 
B = 2Q + W2 - W4, 

c = h 



(4.4) 



^1 — Wb- 



Note that equations of motions (4.3) do not depend on the metric parameter T and on the 
torsion parameter t<i- Moreover, the dependence on W2 and W4 appears only through their 
difference. 

System (4.3) to lowest order becomes 



dv 
~dt 



-L — 2"6 r , 



where e r is the unit vector in the radial direction. Imposing the Newtonian limit it follows 
that (see also formula (23)]) 

C = 1. (4.5) 

Since the Newtonian limit fixes the value of t\, the equations of autoparallels depend only 
on the parameters w\ , . . . , W5 (called frame-dragging torsion parameters in [H] ) . Therefore 
the precession of satellite's orbital elements will depend only on such torsion parameters, 
as it has been found in [9] for gyroscopes. 

Using (4.3) and (3.4) we obtain the following system for the equations of motion: 



( in 



(V + A)xyx + (-Vx 2 + Ay 2 + Bz 2 ) ij + (A-B) yzz 
-(V + A)xyy + (-Ax 2 + Vy 2 - Bz 2 ) x - {A - B) 



j-z : 



(4.6) 



r 2 z+ ejI ^{V + B)z{yx-xy). 



Note that in case of no torsion (i.e. Wi = for any i = 1, . . . , 5) and when Q = —2 system 
(4.6) reduces to the equations of motion found by the Lense-Thirring [151 formula (15)]. 
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5 Computation of orbital elements via perturbation theory 



The system ( |4.6[ ) expressing the motion along autoparallel trajectories can be written in 
the form 



in 



x — ~x -\- F x i 



m 



y + F t 



(5.1) 



■m 



:Z + F Z , 

where (F x , F y , F z ) is the perturbation with respect to the Newton force 
ma 



x — ^5 



F, 



ma 
v = 75" 



(V + A)xyx + (-Vx 2 + Ay 2 + Bz 2 ) ij + (A-B) yzz 
-(T> + A)xyy + (-Ax 2 + Vy 2 - Bz 2 ) x- (A-B) xzz 



(5.2) 



TfiCL 

F z = -^(V + B)z(yx-xy). 

We use the standard coordinates transformation |28| . |29j used in Celestial Mechanics 

x = r (cos u cos — sin u sin Q cos i) , 

y = r (cos u sin Q + sin it cos Q cos i) , 

z = r sin u sin i, 

where i is the orbital inclination, 0, is the longitude of the node, and u is the argument 
of latitude. The vector (F x ,F y ,F z ) can be decomposed in the standard way along three 
mutually orthogonal axes as 

^* y> 



T = d(x/r) F + d(y/r) F + d{z/r) p 



du 



du 



du 



(5.3) 



sinn W = + m^F y + d -^lF z . 

01 01 01 

Here S is the component along the instantaneous radius vector, T is the component per- 
pendicular to the instantaneous radius vector in the direction of motion, and W is the 
component normal to the osculating plane of the orbit (colinear with the angular momen- 



tum vector). Then, substituting (5.2) into (5.3) gives 

J 



S 



T 



W 



jVcosi u, 



J A 

■^Acosi r, 



(5.4) 



J 



sin i (A cos u r — B sin ur u) . 
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Note that in case of no torsion and when Q = — 2 formulae (5.4) reduce to the components 
found by Lense-Thirring (see equations (16) in |15j). 

Let us now recall [28], |29j that, using the method of variation of constants, 

o(l -e 2 ) 



1 + e cos v ' 



where a is the semimajor axis of the satellite's orbit, e is the eccentricity, v is the true 
anomaly, and 



r esmv 



a(l 



2- 2/i 2\l/2 

r v = na (1 — e ) ' , 



n = 2ir/U, U the period of revolution. Following the standard astronomical notation, we 
let oj be the argument of the perigee, and oj = Vt + oj be the longitude of the perigee. 
We also recall the following planetary equations of Lagrange in the Gauss form |29l Ch. 6, 
Sec. 6]: 



" da _ 
~dl ~ 

de 

~dt ~ 
di 

dt ~ 

dn 

~dt ' 

duj 
~dt " 

dL 
dt 



n ( 1 _ e 2)l/2 

(1-e 2 ) 1 / 2 
na 



o(l 



Se sin v + T 



b sin v + 1 e H cos v 



no 2 (l - e 2 )V2 
1 



na 2 (l — e 2 ) 1 / 2 sini 
(l-e 2 ) 1 ^ - 



Wr cos u, 

Wr s'mu, 



(5.5) 



nae 



S cos v + T 1 + 



a(l 



+ 2 sin 2 



2 ~dt' 



Sr + 



na^ 



duj 



l + (l-e 2 )V2 dt 



+ 2(1 



2\l/2 ■ 2 * ^ 

5 'sin — , 

; 2 dt 



where L 



o 



-rn + w is the longitude at epoch, and r is the time of periapsis passage. 



Using the expressions of S, T and W given by (5.4) and integrating the Lagrange planetary 
equations we compute the variations of the orbital elements. According to perturbation 
theory, we regard the orbital elements as approximately constant in the computation of 
such integrals. Since u = v + oj — Q, we can make use of the approximation 



u — v. 



(5.6) 
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Inserting (5.4)-(5.6) into (|5.5l) yields 



da 


2Jecosi (1 + 


~dt ~ 


na 2 (l — 


de 


J cos i sin v 


dt ~ 


na J (l — e*) a i z 


di 


J sin i cos u 


dt 


na 3 (l — e 2 ) 3 / 2 - 


dn 


J sin -u 


~dt 


na 3 (l — e 2 ) 3 / 2 - 


dui 


Jcosi 


~dt 


na 3 e(l — e 2 ) 3 / 2 


dL 
dt ~ 


2 Jcosi _ 
na 3 (l — e 2 ) 



(Av + Vu) , 



e(e + 2 cos w + e cos 2 v)Av + (1 + e cos v) 2 Vu 



(1 + e cos v ) 2 cos uDu — e sin 2 u(2 + e cos v)Av 



+ 2 sin 2 



i dfl 
2 



(1 + e cos v)T>ii + 



l + (l- e 2 )i/ 2 ^ + 2(1 - e ) / Sm 2^ 



(5.7) 

Recalling (5.6), we now integrate (5.7) with respect to v. Therefore we find for the variations 

1 



of the orbital elements: 

2 Je cos i cos v 



da 



, (A + V) ( 1 + ecosv + ^e 2 cos 2 u 



<5e 



J cos i cos t> 



na 



3(! _ e 2)3/2 



(.4 + £>) ( 1 + e cos u + -e 2 cos 2 v ) - A{1 - (?) 



Si 



5Q 



Jsini 



4(A + 2B)e cos v cos 2 u - 4(B + 2A)e cos v 



12na 3 (l - e 2 ) 3 / 2 
+2(0 + 24)e sin v sin(2u) + 30 cos(2n) 

J \ „ 30 . . 

6na 3 (l-e 2 ) 3 / 2 ~ T sm ( 2u ) 



+e 



2(4 - 0) sin u + (4 + 20) sin(2u) cos v - 2(2A + 0) sin u cos 2 it 



SLn 



J 



3 e(l - e 2 ) 3 / 2 



cosz < smu 



V + (4 + £>)ecosu + -(22? - A)e 2 



+-(4 + £>)e 2 cos 2 ?; 



+ (V-A)ev } +2 sin 2 - 5Q, 



2 Jcosi . . e 

-I? (u + esinw) + 



na 3 (l — e 2 ) 



i + (i _ e2) - I7 ,^ + 2(l-e 2 )V 2 s in 2 ^a 
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We note that the contributions of the components S and T to the derivative % are pro- 
portional to Dii and Av, respectively, with the same proportionality constant. Using the 
approximation u c± v it turns out that in the classical Lense-Thirring case, where the torsion 
parameters vanish and —A = T> = Q, there is a cancellation of such contributions in such a 
way that 5a vanishes. Conversely, in presence of torsion, if the eccentricity of the orbit is 
nonzero, the contributions of the radial and of the tangential component of the perturbative 
force differ, so that 5a does not vanish, yielding a periodic perturbation of the semimajor 
axis of the satellite's orbit. 



6 Torsion corrections to the Lense-Thirring effect 



We observe that only periodic terms appear in 5a, 5e and 5i. Secular terms appear in <5f2, 
5uj and 5Lq. Since v = nt + periodic terms in v, the secular contributions to the variations 
of the corresponding orbital elements are: 

J 



(<^)scc 
(5uj) sec 



2a 3 (l - e 2 ) 3 / 2 
J 



Bt, 



a 3(!_ e 2)3/2 
J 



V - A - (B + 2V - 2A) sin 2 



a 3 (l-e 2 ) 



2V + 



l + (l_ e 2)l/2 (l_ e 2)l/2 _ 



V - A - (B + 2V - 2A) sin 2 



-(£ + 4P)sin 2 - j t. 

(6-1) 

In the absence of torsion and when Q = —2, it turns out that (5u>) sec = (6Lo) Bec , as found 
by Lense-Thirring. 



Using (4.4) we rewrite (6.1). For the nodal rate we obtain 

GJ 



(5fl) s 



S(!_ e 2)3/2 



(l + Mi) t, 



and for the longitudinal rate of the perigee 

2QJ 



(5u) 



j3(l _ e 2)3/2 



1 + M2 — 3(1 + H3) sin 



Since u> = £1 + u, for the rate of the argument of the perigee we find 



(5tu) s 



GJ 



The parameters 



a 3 (l — e 



/'1 



1)3/2 



3 + fix + 2fj, 2 -6(1 + ^3) sin 2 



(6.2) 



(6.3) 



(6.4) 



/'2 



M3 



W2 — W4 

2G ' 

2u>i — W3 + U>5 
^2~G ' 

4wi — IV2 — 2u>3 + Wi + 2^5 
-6G 
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measure deviations from GR. Indeed, when there is no torsion we have W{ = for i = 
1, . . . , 5. When, in addition, Q = —2 the metric is the weak field approximation of a Kerr- 
like metric, and [i\ = n% = [13 = and we get the classical Lense-Thirring formulae [T5] . 
We also give the expression for the rate of the longitude at epoch, namely 



(SL ) 



where 



2QJ 



1 



i 3 (l 



l + (l_ e 2)l/2 (l_ e 2)l/2 



1 + M2 - 1 + A*4 



(6.5) 



+ 



1 + Mi + 



3e 2 



1 



l + (l-e 2 ) 1 /2(i_ e 2)i/2 



1 + ^3 1+2(1 + ^*4 



sin 



/i 4 



Note that /ui, ... ,/i4 do not depend on t\,t2-,J r ,'H. 



7 Torsion corrections to the geodetic precession 

The secular perturbations of the orbital elements computed in the previous sections are not 
the only torsion induced perturbations that are expected. Indeed, a further contribution due 
to solar perturbation is present, namely the geodetic precession in presence of torsion. The 
corresponding perturbations of the orbital elements have been computed in the companion 
paper [16] and they depend only on the torsion parameters ij. 

Since we are interested in putting constraints on the frame-dragging torsion parameters 
wi, ... ,W5, there is a relevant difference between the case of GPB gyroscopes considered in 
[9] and the present problem of orbits of satellites. In [9] the average gyroscope precession 
rate is expressed as 




where S'o is the angular momentum of the spinning gyroscope measured by an observer 
comoving with its center of mass, and the vector f2 e g- of the angular precession rate is a 
linear combination of wo (the orbital angular velocity vector of the gyroscope) and 
(the rotational angular velocity vector of the Earth around its axis). In O e g the coefficient 
of wo is a linear combination of the parameters ti, while the coefficient of oje is a linear 
combination of the parameters Wi. Since the GPB satellite has a polar orbit the vectors u)q 
and uje are orthogonal. The contribution to the average precession due to ujo is the geodetic 
precession of the gyroscope, while the contribution due to loe is frame- dragging, both in the 
presence of torsion. Therefore, in the GPB experiment [10], when measuring the projections 
of the average precession rate of a gyroscope on the two corresponding orthogonal directions, 
it turns out that the linear combinations of the ti and of the Wi torsion parameters can be 
constrained separately. 

On the other hand, in the case of orbital motion of satellites, in the presence of torsion the 
geodetic precession and the Lense-Thirring effect are superimposed as it happens in GR, in 
such a way that the precessions of the orbital elements are simultaneously influenced by both 
effects. In |16j it has been found that the contribution of geodetic precession depends on 
a linear combination of the torsion parameters ti, while the contribution of frame-dragging 
computed in the previous sections depends on a linear combination of the parameters Wi. 
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It turns out that the precession of orbital elements (such as the node and the perigee) both 
depend on ti and Wi, in such a way that without a knowledge of the dependence of such 
precessions on U, it is not possible to put constraints on the Wi. The knowledge of the 
dependence on ti corresponds exactly to the knowledge of the geodetic precession of the 
orbital elements in presence of torsion. 

In GR it is known that the geodetic precession is independent of the orbital elements of the 
satellites (and therefore it is the same both for LAGEOS and the Moon). This property is 
used in GR in order to compute an upper bound to the uncertainty in modeling the geodetic 
precession, and in order to show that the result is negligible with respect to the uncertainty 
in the measurement of the Lense-Thirring effect (see [3J, Supplementary Discussion). Such 
a result is important in order to extract the Lense-Thirring effect from LAGEOS data, and 
it is achieved thanks to the precision of the measurement of geodetic precession by means 
of lunar laser ranging (LLR) data [30] . 

In Section [8] we will show that the uncertainty in modeling the geodetic precession can be 
neglected also in presence of spacetime torsion. In particular, the upper bounds on the 



torsion parameters U found in [16] and recalled in the subsequent formula (7.4) will be 
useful in order to obtain such a conclusion. This is important in order to extract the Lense- 
Thirring effect from LAGEOS data also in the presence of torsion, and that will allow us 
to constrain suitable linear combinations of the parameters Wi separately. Hence, in the 
following we briefly need to report the results obtained in |16j . 

The geodetic precession of orbital elements of the satellite in the gravitational field of the 
Earth and the Sun (both supposed to be nonrotating) is computed, in a Sun-centered 
reference system. It is shown that, to the required order of accuracy, the corresponding 
metric is described by a further parameter X = 2(/3— 7), where (3 is the usual PPN parameter, 
and the parametrization of the torsion tensor involves a further parameter £3 (see [16] for 
the details). 

The secular contributions to the precessions of the node and of the perigee due to torsion 
found in [16] are the following: 



(^ SUn )sec=^(C 1 -C 2 ^COSz) t, 

4 p \ n J 

(<5£ Sun ) sec =\ — { Ci + C 2 — [4 - cos i - 5 sin 2 i sin 2 (5 - 0)1 ) t, 
4 o I n i 



(7.1) 



where 



Ci = 1- y + 2X + 3i 2 , C 2 = l + ^ + ^--F-Z + t 2 + 2h. (7.2) 

Here M is the mass of the Sun, vq is the revolution angular velocity of the Earth around 
the Sun, and p is the distance of the Earth from the Sun. 



Differently from the Lense-Thirring effect, the precessions (7.1) depend on the torsion pa- 



rameters t 2 and £3, and are independent of t±\ the parameter t\ is identified using the 



Newtonian limit (4.5). 



We recall that t% and £4 enter the parametrization of torsion at the higher order of accuracy 



required in the computation of precessions (7.1) 



The perturbations (7.1) have to be superimposed to the ones computed in Section 6 



The first term on the right hand sides of the two formulas in (7.1 ) can be interpreted as the 
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geodetic precession effect, when torsion is present [E]: accordingly we set 

(5W C °)sec = (<^ ge °)sec = ~^y^ (7-3) 

In the PPN formalism we have 

Ci = 2 + 4 7 + 3t 2 , C 2 = t 2 + 2(1-/3 + t 3 ). 

Using LLR data and Mercury radar ranging data respectively, the following upper bounds 
are given in [HI Section 13]: 

|t 2 |< 0.0128, |1 - + t 3 \ < 0.0286. (7.4) 

Since for LAGEOS satellites ^ ~ 4.2 x 10~ 4 , we have 

(<m Sun ) scc ~ (^ Se °)sec, (^ Sun ) scc * (^ eo ) sec . 

Taking into account the expression of C±, we have 

(<ms eo )sec = (^ eo ) sec = ^ (l + 2 7 + ht 2 \ t. (7.5) 

This formula yields the rate of geodetic precession around an axis which is normal to the 
ecliptic plane. The projection of this precession rate on the axis of rotation of Earth is 
obtained by multiplying (<5f2 geo ) sec by cos e, where e = 23.5 degrees is the angle between the 
Earth's equatorial plane and the ecliptic plane [31] : this gives the values of the geodetic 
precession in a Earth-centered reference system. 

8 Constraining torsion parameters with LAGEOS 

In this section we describe how the LAGEOS data can be used to extract a limit on the 
torsion parameters. We will assume in the following that all metric parameters take the same 



form as in the PPN formalism, according to (10.2). Recent limits on various components of 



the torsion tensor, obtained in a different torsion model based on the fact that background 
torsion may violate effective local Lorentz invariance, have been obtained in [32] . See also 
[33] . where constraints on possible new spin-coupled interactions using a torsion pendulum 
are described. 

8.1 Constraints from nodes measurement 

Here we discuss how frame dragging torsion parameters can be constrained by the measure- 
ment of a suitable linear combination of the nodal rates of the two LAGEOS satellites. 



Equation (6.2) can be rewritten as 



<«w = a , (1 ! J e2)3/2 ( - \ - ^) t = ««. (") 

where we have defined, similarly to [9] and [16], a multiplicative torsion "bias" relative to 
the GR prediction as 



n = wm = " 2 " = 2 1 1 + 7 + t 1 - (8 - 2) 
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(<5f2)^ = a - i( ^ J e 2^z/2 ^ being the Lense-Thirring precession in GR. We recall that the values 
of such precessions are 31mas/yr and 31.5mas/yr for LAGEOS and LAGEOS II, respec- 
tively, where mas/yr denotes milli-arcseconds per year. 

Let us now consider the contribution of the geodetic precession to the nodal rate. We write 
the secular contribution to the nodal rate, in a Earth-centered reference system, in the form 

(e^ e °) S ecCOSe= (^eo )s GR cose 6 gec ; ^ 

where 6q CO depends on t<i- Precisely, taking into account that (5r2 geo )^ = 3J j^ / ° 1 and using 



(7.5), we have 

6^ = i(l + 2 7 ) + |. (8.4) 

Moreover, the following numerical constraints are set on PPN parameters 7 and ct\ by 
Cassini tracking [33] and LLR data [3], respectively: 

7- 1 = (2.1 ±2.3) x 10~ 5 , |ai|<10" 4 . (8.5) 



From ( 
using 



8.5) it follows that the term i(l±2 7 ) differs from 1 by a few part in 10 5 . Therefore, 



7.4), (8.4) and (8.5) we get 



l&n°-l| - if I < 0.0064. (8.6) 

The measurement of the Lense-Thirring effect in [3J, [35] is based on the following linear 
combination of the total nodal rates of the two LAGEOS satellites: 

sn\ ot + K8n\?, (8.7) 

where the subscripts I and II denote LAGEOS and LAGEOS II, respectively. Here the 
total nodal rate 5Q tot of a LAGEOS satellite denotes the nodal rate due to all kinds of 
perturbations, both gravitational and nongravitational. The coefficient k = 0.545 is cho- 



sen to make the linear combination (8.7) independent of any contribution of the Earth's 
quadrupole moment J2, which describes the Earth's oblateness. 

In [3] the residual (observed minus calculated) nodal rates A(50i), A(<5f2n) of the LAGEOS 
satellites are obtained analyzing nearly eleven years of laser ranging data. The residuals 
are then combined according to the linear combination A(5Cli) + kA(5CIu), analogue to 



(8.7). The Lense-Thirring effect is set equal to zero in the calculated nodal rates. The 
linear combination of the residuals, after removal of the main periodic signals, is fitted with 
a secular trend which corresponds to 99% of the theoretical Lense-Thirring prediction of 
GR (see [3J, [33] for the details): 

(<5Ql) s G * + K(5n n )?£ = 48.2 mas/yr. 

The total uncertainty of the measurement is ±5% of the value predicted by GR [3J, [35] . 
[36] . This uncertainty is a total error budget that includes all estimated systematic errors 
due to gravitational and non-gravitational perturbations, and stochastic errors. Such a 
result is quoted as a 1 — a level estimate in [37] , [3H] , though an explicit indication of this 
fact is missing in [4J. Eventually, the authors allow for a total ±10% uncertainty to include 
underestimated and unmodelled error sources [4]. In the following we assume a value of 
±10% for the uncertainty of the measurement. 
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Using the upper bound (8.6), the uncertainty in modeling geodetic precession in the presence 
of torsion is 



geoxGR 
/sec 



cose < 0.0064 



27.2 r 



48.2 



GR 

sec 



where [(<5f2j eo )^ + ft^f^ )^] cose = 27.2 mas/yr is the contribution from geodetic pre- 
cession predicted by GR for LAGEOS satellites. Compared to the ±10% uncertainty in 
the measurement of the Lense-Thirring effect, the uncertainty in modeling geodetic preces- 
sion can be neglected (as in [3j, [35]) even in the presence of spacetime torsion. This is a 
consequence of the torsion limits set with the Moon and Mercury in [16] . 
Then we can apply the results of [3], |35j to our computations with torsion, and we obtain 



(<^i) sec + K(6n n ) sec - o.99 [{SQi)^ + «(<yn n )S 



< o.io [(<yni)™ + K(sn n 



)G R1 

/sec J j 



where (50i) sec and (tffln)sec are given by (8.1). Since the torsion bias 6q does not depend 
on the orbital elements of the satellite, we have 



(£Ol) B ec + ft(£On) S ec 



bn. 



Hence, using (8.2), we can constrain a linear combination of the frame-dragging torsion 
parameters W2, W4, setting the limit 

1 



|6n — 0.99] 



7- 1 + 



V)2 — W4 



+ 0.01 



< 0.10, 



which is shown graphically in Figure 1, together with the other constraints on 7 and ot\ [3]. 
Taking into account the numerical constraints (8.5) the limit on torsion parameters from 
LAGEOS becomes 

W2 — 



+ 0.02 



< 0.20 



which implies 



The constraint 



0.36 < w 2 - w A < 0.44. 



(8.9) 



3.9) on the torsion parameters depends on the quantitative assessment of 
the uncertainty of the measurement of the Lense-Thirring effect. However, the value 5-10% 
of the uncertainty reported in [1] has been criticized by several authors. For example Iorio 
argues in |37J that the uncertainty might be 15-45%. The previous computations show that 
the upper bound on the quantity 



W2 — W4 



+ 0.01 



is given by the uncertainty of the measurement, so that one can find the constraint on the 
linear combination of the torsion parameters u>2 , W4 corresponding to a different value of 
the uncertainty. For instance, if the value of the uncertainty of the measurement is ±50%, 
the constraint on torsion parameters becomes 

-1.96 < w 2 -wa < 2.04. 

One of the goals of the LAGEOS, LAGEOS II, LARES^ three-satellite experiment, together 
with improved Earth's gravity field models of GRACE (Gravity Recovery And Climate 



J LAser RElativity Satellite, a geodynamics mission of the Italian Space Agency (ASI) to be launched. 
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Experiment) is to improve the experimental accuracy on the orbital Lense-Thirring effect 
to "a few percent" [35] . 



We observe that, using (8.8) the uncertainty in modeling geodetic precession in presence of 



torsion amounts to about 0.4% of the Lense-Thirring effect, which is still a small contribu- 
tion to a total root-square-sum error of a few percent. Note that an improved determination 
of the geodetic precession has been recently achieved by GPB [TDj which, unlike LAGEOS, 
is designed to separate the frame-dragging and geodetic precessions by measuring two dif- 
ferent, orthogonal precessions of its gyroscopes. 




Figure 1: constraints on PPN parameters (7, a±) and on frame-dragging torsion parameters 
(1V2, W4) from solar system tests. The grey area is the region excluded by lunar laser ranging 
and Cassini tracking. The LAGEOS nodes measurement of the Lense-Thirring effect [5J, |35j 
excludes values of (W2 — w±) /2 outside the hatched region. General Relativity corresponds 
to 7 = 1, a\ = and all torsion parameters = (black dot). 



In the case of GPB, the torsion bias for the precession of a gyroscope is (9] 

Q W\ + W2 — U>3 — 2ui 4 + w§ 
~ 2 2 ' 

This formula (the analogue of the right hand side of equation ( |8.2| )) involves a linear com- 
bination of all frame-dragging torsion parameters. Such a linear combination can be con- 
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strained from GPB data. Since LAGEOS and GPB are sensitive to different linear combi- 
nations, together they can put more stringent torsion limits. 

After taking into account the contribution of the geodetic precession, the combined con- 
straints from gyroscope and orbital Lense-Thirring experiments are effective probes to search 
for the experimental signatures of spacetime torsion. In this sense, LAGEOS and GPB are 
to be considered complementary frame-dragging and, at the same time, torsion experiments, 
with the notable difference that GPB measures also the geodetic precession. 



8.2 Constraints from nodes and perigee measurement 



In this section we discuss how frame dragging torsion parameters can be constrained by the 
measurement of a linear combination of the nodal rates of LAGEOS and LAGEOS II and 
the perigee rate of LAGEOS II. 

Similarly to the previous section, we define a multiplicative torsion "bias" relative to the 



GR prediction also for the rate of the argument of the perigee (6.4): 



(5u>) 



GR 

sec 



6 cos i 



3 + Mi + 2^2 - 6(1 + /Z3) sin 2 



(<5w)sec = — °2^ Z /2 ~k being the Lense-Thirring precession in GR: we recall that the value 



of this precession is — 57mas/yr for LAGEOS II. In the following, the torsion bias b u 
referred to LAGEOS II. 

Using the values of fj,±, fj>2 and ^3 given in Section [6] we find 



is 



"2 + 



12 ' 



UO) 



The measurement of the Lense-Thirring effect in |39| is based on the following linear com- 
bination of the residuals of the nodes of LAGEOS and LAGEOS II and of the perigee of 
LAGEOS II: 

A{5Qi) + ciA(<5fin) + c 2 A(5uj n ), (8.11) 
where the coefficients c\ = 0.295 and C2 = —0.35 are chosen to make the linear combination 



(8.11) independent of the first two even zonal harmonic coefficients J 2 and J4, and of their 



uncertainties. 

In [39] the residuals are obtained analyzing four years of laser ranging data, and then 



combined according to the linear combination (8.11). The Lense-Thirring effect is set equal 



to zero in the calculated rates of the nodes and of the perigee. The linear combination of 
the residuals, after removal of the main periodic signals and of small observed inclination 
residuals, is fitted with a secular trend which corresponds to 1.1 times the theoretical Lense- 
Thirring prediction of GR (see [39] for the details): 



(50 T ) S G * + Cl (£fi„) s G * + c 2 (fon) s G * = 60.2 mas/yr. 

The total uncertainty of the measurement found in [39j is ±20% of the value predicted 
by GR. This uncertainty is a total error budget that includes all the estimated systematic 
errors due to gravitational and non-gravitational perturbations. Such a result is quoted as 
a 1 — a level estimate in [37] , though an explicit indication of this fact is missing in [39J . 
The contribution to the uncertainty of the measurement due to nongravitational perturba- 
tions, mainly thermal perturbative effects, on the perigee of LAGEOS II, amounts to 13% 
of the value predicted by GR. In |40j such an estimate is confirmed, however the author, 
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when considering more pessimistic assumptions on some thermal effects, estimates that the 
contribution of nongravitational perturbations to the total uncertainty does not exceed the 
28% of the GR value. Here we will follow this more conservative estimate. Inserting this 
value in the estimate of the total uncertainty computed in [39] yields a total root-square-sum 
error of 32% of the GR value. 

For reasons similar to the ones discussed in the previous section, we are allowed to neglect 
the uncertainty in modeling the geodetic precession in presence of torsion. Then we can 
apply the results of [39] to our computations with torsion, and we obtain 

(^i) scc + Cl (<5Q„) sec + c 2 (^ii) sec - 1.1 [(<ft2i)Sc + ci(<5fin)™ + c 2 ((^n) s G *] 



< 



A direct computation gives 



0.32 [(«2i)<^ + ci(^n) s G j? + c^n)^] 



[l-K)ba + Kb u -l.l\ <0.32, 



1.12) 



where 



K 



c 2 (5uju) 



GR 

sec 



0.33. 



Inserting in (|8.12|) the expressions of ba and b^ given in (|8.2|), (|S. 10[) and taking into account 
that Q 



-2 by formula (8.5), we obtain 

U>4 



-0.22 < 



U>2 



+ K 



2lU\ + U>2 — W3 — W4 + W5 

6 



< 0.42. 



Using the value of K we finally deduce 



0.22 < O.lliui - 0.20^2 - 0.06^3 + 0.20w 4 + 0.06u> 5 < 0.42, 



5.13) 



which is shown graphically in Figure 2, together with the other constraints on 7 and a\ [3J. 



The constraint (8.13) on the linear combination of the frame-dragging parameters is rather 



weak, due to the uncertainty on the nongravitational perturbations. Notice that the coef- 
ficients in front of W3 and W5 are of an order of magnitude smaller than the coefficients of 
the other parameters, so that the constraint on u>3 and W5 is even looser. 
Thermal thrusts (TTs) are the main source of non-gravitational perturbations [3Q]. One of 
the main drivers of LAGEOS TTs is the thermal relaxation time tccr of its fused silica 
cube corner retroreflectors [H] , which has been characterized in laboratory-simulated space 
conditions at the INFN-LNF Satellite/lunar laser ranging Characterization Facility (SCF) 
|42j . |43j . |44j . The measurements of LAGEOS tccr in a variety of thermal conditions 
provide the basis for possibly reducing the uncertainty on the thermal perturbative effects. 
As a consequence, the constraint 



.13) could be improved. 



The constraint (8.13) on the torsion parameters depends on the quantitative assessment of 



the uncertainty of the measurement of the Lense-Thirring effect. Again, the value ±20% 
of the uncertainty reported in [39] has been criticized by various authors. For example 
Ries, Eanes and Tapley argue in [45J that the uncertainty is at best in the 50-100% range. 
The uncertainty of the measurement yields the upper bound on the right-hand side of the 



estimate (8.12). Hence, one can find the constraint on the linear combination of the torsion 



parameters Wi corresponding to a different value of the uncertainty as it has been discussed 
in Section (HO 
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Figure 2: constraints on PPN parameters (7, a\) and on frame-dragging torsion parameters 
(wi, u>2, W3, W4, W5) from solar system tests. The grey area is the region excluded by lunar 
laser ranging and Cassini tracking. The LAGEOS nodes and perigee measurement of the 
Lense-Thir ring effect [39], [30] excludes values of 0.1 lwi — 0.20^2 — O.O6W3 + 0.20^4 + 0.06^5 
outside the hatched region. General Relativity corresponds to 7 = 1, a± = and all torsion 
parameters = (black dot). 



We recall that in [9] an upper bound on the combination \w\ + W2 — W3 — 2w^ + 105 1 is 
given. This constrains the torsion parameters within two parallel hyperplanes in a five- 
dimensional space. If we couple this bound with our two estimates (8.9) and (8.13), we 
obtain that wi, ... ,11)5 are constrained to lye in a five-dimensional set, which is unbounded 
only along two directions. Hence, coupling GPB with SLR measurements significantly 
reduces the degrees of freedom on the frame-dragging parameters. 

We conclude this section by observing that the recently approved JUNO mission to Jupiter 
|46| will make it possible, in principle, to attempt a measurement of the Lense-Thirring 
effect through the JUNO's node, which would be displaced by about 570 metres over the 
mission duration of one year [37]. Hence, such a mission yields an opportunity for a possible 
improvement of the costraints on torsion parameters. 
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9 Conclusions 



We have applied the framework recently developed in [9j for GR with torsion, to the compu- 
tation of the slow orbital motion of a satellite in the field generated by the Earth. Starting 
from the autoparallel trajectories, we computed the corrections to the classical orbital Lense- 
Thirring effect in the presence of torsion. By using perturbation theory, we have found the 
explicit dependence of the secular variations of the longitudes of the node and of the perigee 
on the frame-dragging torsion parameters. The LAGEOS nodes measurements [3], |35] and 
the LAGEOS nodes and perigee measurements [39], [30] of the Lense-Thirring effect can be 
used to place constraints on torsion parameters, which are different and complementary to 
those set by GPB. 



10 Appendix 

Under spherical axisymmetry assumptions, the metric tensor can be parametrized to 
first order as follows [31: 



ds 2 



i+n- 



. m 



dt 2 + 



dr 2 + r 2 (d9 2 + sin 2 9 d<p 2 ) + 2Q- sin 2 9 dtdcp, (10.1) 

r 



where H,J-,G are three dimensionless parameters that can be immediately related to the 
Parametrized Post Newtonian (PPN) parameters: 



n 



■F = 2 7 , 



1 + 7- 



cti 



(10.2) 



Here we follow the notation of the paper [9j, instead of the PPN notation. This will be 
useful in Section 

The nonvanishing components of the torsion tensor are: 



q t 



s, 



rB 



q t 



H4> 



'tr 



J t0 



hi 

2T 2 ' 



t2 



m 
2^2 



J ■ n 

W2 — sin 9 cos I 

2r 

^3^2 sin 2 ^, 



(10.3) 



W4— r sint^ cos( 
2r A 



1V5 



J 

2r 3 
J 

2^' 



J cos 9 



-W4- 



' 2r 3 sin 9 

The expression of the nonvanishing components of the connection approximated to first 
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order in e m = m/r, ej = J/(mr) and e m ej = J/r is the following [9]: 

r*tr = ^~ (2*i - %) e m , 



r r 0t = -(Q - {wi +w 3 + w 5 )) sin 2 8e m ej, 



r* - 

L rt - 2r m ' 



T 4 ^ = ^(3G + (wi - w 3 - w 5 ))sin 2 6 e m ej, 

r*^ r = — (3^ — (wi + w 3 + it's ) ) sin 2 9 e m ej , 

= ^w 2 r sin 9 cos 8 e m e j, 

T 

1 rr — 2 r m ' 

= -r sin 2 6* + ^(J 7 + £2) sin 2 9e m , 

^ r t<t> = \{Q ~ ( w i ~ w 3 + w 5 )) sin 2 9e m ej, 



T e t(j) = --^(20 + (w2-2w A ))sm9cos9e m ej, 

T% = --^(20 + w 2 )sm9cos9e m ej, 

r e _ r _ 1 

r e _ v 4> _ 1 ly. , 

T 8 ^ = - sin 6> cos 0, 

r^ r = -^(^- {wx - w 3 + w 5 ))e m ej, 

rtt = -^(G - (wi -w 3 -w 5 ))e m ej, 

^te = 2 ^( 2 G + (w2-2w 4 )) ( ^e m ej, 

2r sm 



cos 6 1 



sini 
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